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ABSTRACT 

The large-scale structure (LSS) in the Universe comprises a complicated filamentary 
network of matter. We study this network using a high-resolution simulation of struc- 
ture formation in a A Cold Dark Matter cosmology. We investigate the distribution 
of matter between neighbouring large haloes whose masses are comparable to mas- 
sive clusters of galaxies. We identify a total of 228 filaments between neighbouring 
clusters. Roughly half of the filaments are either warped or lie off the cluster-cluster 
axis. We find that straight filaments, on the average, are shorter than warped ones. 
Close cluster pairs with separation of 5/i~^Mpc or less are always connected by a 
filament. At separations between 15 h~^ Mpc and 20 h~^ Mpc about a third of cluster 
pairs are connected by a filament. On average, more massive clusters are connected to 
more filaments than less massive ones. This finding indicates that the most massive 
clusters form at the intersections of the filamentary backbone of LSS. For straight 
filaments, we compute mass profiles. Radial profiles show a fairly well-defined radius, 
Ts, beyond which the profiles follow an r~^ power law fairly closely. For the majority 
of filaments, Vg lies between 1.5ft.~^Mpc and 2.0 /i~^ Mpc. The enclosed overdensity 
inside r^ varies from a few times up to 25 times the mean density, independent of 
the length of the filament. Along the filaments' axes, material is not distributed uni- 
formly. Towards the clusters, the density rises, indicating the presence of the cluster 
infall regions. Filaments have been suggested to cause possible alignments between 
neighbouring clusters. Looking at the nearest neighbour for each cluster, we find that, 
up to a separation of about 15 h~^ Mpc, there is a filament present that could account 
for alignment. In addition, we also find some sheet-like connections between clusters. 
In roughly a fifth of all cluster-cluster connections where we could not identify a fila- 
ment or sheet, projection effects lead to filamentary structures in the projected mass 
distribution. 

Key words: cosmology: theory, methods: N-body simulations, dark matter, large- 
scale structure of Universe 



1 INTRODUCTION 

Large galaxy redshift surveys sucli as the Sloan Digital Sky 
Survey (York et al. 2000) and tlie 2dFGRS (CoUess et al. 
2001) and N-bo dy simulations of cosmic structure forma- 
tion (for example rjenkinsl(ll998l ). IWambsganss et al.lll2004l) ') 
show a complicated network of matter. In redshift surveys, 
the galaxies line up preferentially around roundish, almost 
empty regions - so-called voids. Clusters of galaxies have 
very prominent positions in the network: They lie at the in- 
tersections of filaments. In N-body simulations, this trend 
is even more pronounced. The network of haloes is clearly 
followed by the more diffusely distributed component of the 
dark matter that, at z = 0, has not collapsed into haloes. 

Describing this network is no easy task. There are 
a wide variety of statistical and topological tools to 
compare observations and theoretical models. Amongst 



these are the t wo-point (and higher-order) cor r elatio n 
function(s) (e.g. iPeebleJ Il980l: iPeebles fc Gro th 197^, 
minimal spanning trees (see e.g ^_|Barro^e^l- .19851 : 
iKrzewina fc Saslawl [199^: iBhavsar fc Splinted Il99d) . the 
genus statistics IIGott et alj ^86|) , M inkow ski functionals 
(that inclu de the genu s: iMec ke et al. |[l994i), shape statis- 
tics (see e. g. [Ba bul fc Starkma n 1992; Luo fc VishniaJ 19951 : 
ILuo et aD Il99d) . and shapefinder statistics derived from 
Minkowski functionals llSahni et al.lll998l: ISheth et alj|2003l : 
IShandari n et alJl2004|) . These tools have all been applied to 
galaxy catalogues and N-body simulations and have been 
very useful in comparing observations and theory. 

These tools have been somewhat less helpful in describ- 
ing the pattern of large-scale structure as far as filaments 
and sheets/walls are concerned. Minkowski functionals and 
especially their shapefinder cousins have been used to try to 
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measure how thick or long filaments are. By construction, for 
simple toy models like ideal filaments (cylinders) or sheets 
(planes), shape statistics give very simple answers. Because 
of the complexity of the network of matter in surveys or 
simulations the answers for these cases are usually more dif- 
ficult to interpret than in the toy model simulations. Re- 
cently, |BliMadwaj_et_alJ ll2004f) investigated filaments in the 
Las Campanas Redshift Survey using shapefinders. Adding 
statistical tools to their machinery, they looked at the max- 
imum length scale at which filaments in that survey are sta- 
tistically significant. They found that scale to be 50 /i~^ Mpc 
to 70 h~^ Mp c. Another very pro mising approach has been 
suggested by IStoica et al.l i2004f) . They introduce another 
method to actually detect filaments. Their method is based 
on an algorithm that has been successfully used to detect 
road networks, and they showed that it works very well for 
2D mock data sets. 

In this work, we will restrict our focus mainly to fil- 
ament s. In slices th rough N-body simulations (see for ex- 
ample |jenking|l993) filaments appear to be very common. 
There is some variety in their sizes and appearance. The 
most massive haloes are usually connected by very promi- 
nent filaments (see Figure^ . However, some of the filaments 
could be sheets and only appear to be filamentary due to 
projection effects. Less dense and less massive filaments can 
be found in less dense regions - tho se filaments r esemble fine 
perl necklaces. In an earlier work iColberg et al. 1999) . we 
showed that the formation of clusters is intimately linked 
with the cosmic neighbourhood. Infall of matter into the 
region where a cluster is forming happens from a few pre- 
ferred directions. These agree with the locations of filaments 
at z = 0. 

Here, we will investigate the properties of filaments in 
more detail. In particular, we will address the following set of 
questions: On average, how many filaments intersect at the 
location of a cluster? Does the number of filaments depend 
on the mass of the cluster? What is the typical density of 
a filament? How long are these filaments typically? For a 
pair of clusters at some separation what is the likelihood of 
finding a filament between them? 

This paper is organized as follows: After a brief discus- 
sion of earlier theoretical work ||5J, in Section |3 we discuss 
the status of observations to locate filaments. Afterwards 
(Section 0, we briefly describe the simulation that we use 
(Section 14.111 and then the procedure (Section 14.311 to look 
at cluster-cluster connections. We present a classification 
of those connections in Section 14.41 Section 14.51 describes 
properties of filaments and of clusters connected to such fil- 
aments. In Section fl.Sl we also discuss what our results might 
mean for alignments of clusters. We conclude the paper with 
a summary (Section (KJ. 



2 THEORETICAL WORK 

Anisotropic collapse of matter in gravitational instability 
scenarios has been known to lead to the formation of sheets 
and fila ments since the seminal work b y jZel'dovich ( 1970 ) 
(also see IShandarin fc Zeldovichlll989^ . ]ic"kel Jl973l) looked 
at the effect using homo geneous ellip soidal models (also se e 
the work by IWhVte fc S ilk 1979; Eis enstein fc Loeblll99El) . 
iBond et alJ il996l) then emphasized the role of tidal fields 



and of anisotropic collapse in the formation of LSS. They 
connected tidal shear directly to the locations of filamentary 
structures and of peaks. 

That same year, Ivan de Weveaert fc Bertschingen 
(199^ investigated the typical morphology of a config- 
uration of two clusters with material in between to find 
that the primordial shear constraint naturally evolves into 
a config uration of two clusters that are connected by a 
filament. IVan de Wevgaerd 112002^ contains a very detailed 
summary of these theoretical efforts. 

There have also been theoretical studies of gas in 
filaments, the so-called warm-hot intergalactic medium 
(WHIM), using simulations. The WHIM contains a sig- 
nificant fraction of all baryons in the present-day uni- 
verse (about 30% to 40%), which makes observing fila- 
ments through the signature of their baryons quite inter- 
esting. For the most recen t studies see Dav e et al. (200l|), 
iKravtsov et alj l|2003), and lFurlanetto et alJ (|2003) and ref- 
erences therein. 



3 OBSERVATIONAL STATUS 

The questions we posed in the introduction are of par- 
ticular interest given observational efforts to find fila- 
ment s by lookin g at t he sp ace in the vic i nities of clus- 
ters. iGrav et al.l (|20o3) and iDietrich etaP \2m4\ looked 
at cluster pairs A901/A902 and A222/A223, respectively. 
Using weak lensing they both conclude that the two clus- 
ters i n their respective surveys are connected by a fila- 
ment. |pimbbier^[DrinkwateJ(|20ol) used clusters A1079 
and A1084 in a pilot study to look for filaments, reporting 
a "filament detection at a 7.5(7 level." iTittlev fc HenriksenI 
( 2001) used x-ray data to detect a filament between clusters 
A3391 and A3395. 

The most relevant detection of a filame nt that was 
not fo und by looking at clusters directly is bv IScharf et alJ 
l|2000|). They found a "5o" significance half-degree filamen- 
tary structure" , present both in x-r ay and optica l data . At a 
likely redshift of around >z = 0.3. IScharf et al.l ll2000l) give 
the length of the structure as > 12 h^^ Mpc. 

Superclusters are very likely locations of fila ments or 
even sheets. For example, IConnoUv et alj (|l99w) found a 
large structure at a redshift oi z = 0.54 with an overdensity 
(in galaxies) of about four that includes three X-ray clusters. 
The galaxies in this sample form "a linear structure pass- 
ing from the Southwest of the survey field throu gh to the 
Northeast." More recentlv. iBregman et alJ l|2004fl looked at 
filaments in superclusters through UV absorption line prop- 
erties of three AGNs projected behind possible filaments in 
superclusters. They conclude that their results are consis- 
tent with the presence of filaments. ^_^ 

At somewhat larger redshifts, iGal fc LubinI (|200J) in- 
vestigated two clusters at redshifts of a « 0.9 which appear 
to be connected by a large structure. lEbeling et al.l (|2004r 
detected a structure of galaxies extending out from the clus- 
ter MACS J0717.5-I-3745, that is located at a redshift of 
z = 0.55, with a length o f 4 h^f^ Mpc. 

iKaastra et alJ l|2003l ) observed a sample of 14 clusters, 
looking for soft X-ray excess emission. For five of their clus- 
ters they find "a significant soft excess," which they at- 
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tribute to emission from intercluster filaments of the WHIM 
in the vicinity of these clusters. 

All these works are very exciting and indicate that there 
will be many more such projects in the very near future. We 
thus feel that answering the questions about what one can 
expect to find is all the more relevant. 



4 THE INTER-CLUSTER NETWORK OF 
MATTER 

4.1 The simulation 

For this work, we use the ACDM simulation introduced in 



iKauffmann et al.l fl999). The simulation parameters corre- 
spond very closely to what has recently become the standard 
cosmology with 30% of the critical density contributed by 
Cold Dark Matter and the remaining 70% by a cosmological 
constant (or "Dark Energy"). A Hubble constant of /i = 0.7^ 
is used and the model is cluster normalized to ag — 0.9. With 
256"^ particles in a box of size (141.3 h~^ Mpc)^ the simula- 
tion volume is large enough to study large-scale structure 
at high resolution'^. 

4.2 Visual impression 

Figure shows a slice of thickness 10ft~^Mpc through 
the simulation volume. The dark matter distribution was 
smoothed adaptively, and the resulting density field is shown 
using a logarithmic colour scale. The slice shows the network 
of matter, which is quite familiar from simulation work and 
from galaxy redshift surveys. 

We have marked the location of the filament that is 
shown in Figure |3 with a box. In order to emphazise the 
role of the clusters we have included the clusters inside the 
box (Figure |3 shows just the material inside the filaments, 
excluding the clusters). The width of the box - its dimension 
perpendicular to the cluster-cluster axis - is 15^~^Mpc, 
which is the actual diameter of the cylinders used to cut out 
filaments. Figure Ogives a good impression of the extent of 
the longer filaments in our sample'^. 

4.3 The procedure 

First, we extract clusters from the halo catalogue that we 
downloaded from the website of the Max-Planck-Institut 
fiir Astrophysik"*. We consider all haloes more massive than 
lO^'' M0. This provides a sample of 170 clusters whose num- 
ber density is roughly com parable to that of R=0 Abell clus- 
ters JPostman et alJll992ll . 

Since we are interested in the inter-cluster network of 
matter, we examine the twelve nearest neighbours of each 



^ Throughout this work, we will express the Hubble constant in 
units of -Ho = lOOh km/sec/Mpc. 

^ The simulation data and halo catalogue can be down- 
loaded from |http : //www ■ mpa-garching . mpg ■ de/galf orm/virgo/ 1 
hrs/index . shtml. 

^ Note that the large halo that is close to the cluster at the top of 
the box has a mass below our cluster threshold mass. See Section 
1131 for details. 



|http : //www . mpa-garching ■ mpg . de/galf orm/vlrgo/lndex ■ shtml | 



cluster. For each pair of clusters we extract the dark matter 
between them as follows. The cluster centers define an axis. 
We extract all dark matter particles in a cylinder of radius 
7.5/i~^Mpc around the axis. Furthermore, we only work 
with those particles that, when projected onto the axis, lie 
outside both clusters' r2oo- That way, we do not extract 
matter that belongs to either cluster. However, matter that 
lies in the cluster-infall regions is included in the samples. 
The value of 7.5 h~^ Mpc is empirical (as is the number of 
twelve cluster neighbours). We found that going further out 
does not add any extra information whereas using smaller 
radii sometimes tended to cut warped filaments in half. 

There is one final requirement for how neighbouring 
clusters were picked. Because we are interested in the inter- 
cluster network of matter we want to avoid finding one clus- 
ter lying directly between two other clusters. We thus do not 
allow another cluster to lie inside the innermost 5h~^ Mpc 
from the cluster-cluster axis. As before, this value is empir- 
ical. 

We want to emphasize that most certainly there will be 
many configurations where a number of clusters are lined up 
on a long filament. We exclude these cases for the simple and 
only reason that we want to study the distribution of matter 
between two neighbouring clusters. In a follow-up study, we 
will come back to looking at the alignment of clusters on 
larger scales. 



4.4 Classifying inter— cluster connections 

Conventional wisdom says that if one picks two neighbouring 
clusters, they are connected by a filament, they lie in a sheet, 
or there is a void between them. In reality, the number of 
possibilities is slightly more complicated. The following list 
of the configuration of matter between neighbouring clus- 
ters is based on the visual inspection of each of the cluster- 
cluster connections. We tried to automate the process by 
using density measures but there are too many cases with 
deviations from the simple configurations mentioned above. 

Filaments: 19% of the 1207 cluster-cluster connections^ 
contain a filament. We found three different possible config- 
urations: 

straight: 38% of all filaments are straight and on center 
with respect to the cluster-cluster axis. That is, the clusters 
lie on the axis of the filament. See Figure |21 

off center: Another 9% of all filaments are also fairly 
straight filaments but their central axes do not align with 
the axis that connects the cluster centers. 

warped/irregular: 53% of all filaments are not straight 
but are either warped or consist of multiple parts. Warped 
filaments sometimes indicate the presence of another cluster 
that lies just outside the 5 h~^ Mpc exclusion zone. When we 
looked at the matter distribution of warped filaments going 
beyond the cylinder used for the analysis we found many 
cases where nearby mass concentrations must have tidally 
interacted with the filaments. See Figure |3 



^ Usually, if cluster rii has cluster n2 in its list of neighbours, 
rii also shows up in n2's neighbour list. However, this is not al- 
ways the case. Thus, the list of connections deviates from the 
completely symmectric case 170 X 12/2 = 1020. 
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Figure 1. A slice of tliickness 10/i~^Mpc tlirougli the GIF simulation. The dark matter was smoothed adaptively, and the resulting 
density field is shown using a logarithmic colour scale. We have marked the location of the filament that is shown in Figure |21 with a 
box. In order to emphazise the role of the clusters we have included the clusters inside the box. The width of the box - its dimension 
perpendicular to the cluster-cluster axis — is 15/i~^Mpc, which is the actual diameter of the cylinders used to cut out filaments. - 
Download fromjhttp://lahmu.phyast.pitt.edu/~colberg/GIF_L0_42_z.jpg| 
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Figure 2. Two orthogonal projections of the dark matter between two of the clusters in the GIF simulation. The plots show the projected 
overdensities, smoothed with a Gaussian of radius 0.5 /i~^ Mpc. The contour levels show overdensities ranging from 0.0 (mean density, 
black) to 19.0 (white). The y-axis cuts through both cluster centers. The region shown here excludes the clusters themselves. Matter 
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Figure 3. Same as Figure |21 for a moderately warped filament. Please note that this cluster-cluster connection is longer than the one 
shown in Figure 1^ 
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Figure 5. Same as Figure |3- again with different length scale - for matter that appears to be filamentary from one direction (left 
panel). Viewed from another direction — on the right-hand side from a direction perpendicular to the one on the left-hand side there is 
no visible filament. The filament on the left-hand side is merely caused by a projection effect. 



Sheets: In about 2% of all cases, we found a sheet-like 
configuration between cluster pairs. Sheets display a narrow 
extent when viewed edge on and a broad uniform distribu- 
tion of matter otherwise. See Figure 0] 

The Rest: The remaining 79% of cluster-cluster pairs do 
not fall into the aforementioned categories. However, there 
are some interesting cases left. In 3% of all cluster pairs we 
found a large amount of matter between them that does not 
look either filamentary or sheet-like. Instead, viewed from 
any angle matter fills the space between the clusters almost 
uniformly. In 19% of all cases that were not classified as a 



filament, sheet or crowded field projection effects lead to the 
appearance of filaments. When viewed from one angle there 
seems to be a filament, whereas viewed from another angle 
there is none (see Figure (SJ. There is the possibility of a 
sheet being mis-classified as a projection effect or vice versa 
- especially since the classification is done by eye. However, 
since we focus our work on studying filaments this uncer- 
tainty does not affect the bulk of our results. Projection ef- 
fects and, to a lesser extent, sheets have to be kept in mind 
when looking at actual measurements of the distribution of 
matter between clusters by means of gravitational lensing as 
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Figure 6. Fractional abundance of filaments: The plot shows the 
fractions of cluster-cluster connections that contain a filament - 
straight or warped — as a function of length of those connections. 
Close cluster pairs are always connected by a filament. 



Figure 8. Enclosed overdensity profile of a straight filament. The 
bold vertical line shows the radius at which the profile starts to 
follow an r~^ power law. 
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Figure 9. Distribution of scale radii r^ of straight filaments. 



Figure 7. Length distributions of filaments: The plot shows the 
fraction of filaments (solid line) and warped filaments (dotted 
line) that have a length I. Warped filaments tend to be longer 
than straight filaments. 



results from these very different configurations will appear 
the same. 



4.5 Properties of inter— cluster filaments 

Our sample of straight and warped filaments is large enough 
to allow more detailed studies of their properties. It is im- 
portant to bear in mind the restrictions imposed on the anal- 
ysis by how we selected the filaments. By filament we mean 
a filament that connects two clusters. The filament could be 
a segment of a larger filament that connects more than two 
clusters. We do not attempt to study large filaments in more 
detail here. Whenever we are talking about filaments in the 
following we mean filaments between pairs of clusters. 
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Figure 10. For straight filaments, this plot shows the enclosed 
overdensity at scale radius rs versus the length of the filament. 
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Figure 11. Longitudinal density profile of straight filaments, av- 
eraged over straight filaments that are longer than 5/i~^Mpc. 
Shown is the enclosed overdensity as a function of the positions 
along the cluster-cluster axis for all material that is contained 
within 2 h~^ Mpc from the axis. 



4.5.1 General properties of cluster-cluster filaments 

The first questions to ask about filaments are how long and 
how common they are. As we noted before, we looked at 
1207 regions between neigiibouring clusters of which 19% 
contained a filament. We selected the pairs of clusters by 
picking the twelve nearest neighbours of a given cluster. 
Given the fact that in the large-scale structure matter sur- 
rounds large voids, the relatively small number of filaments 
does not come as a surprise. If we had required that the 
cluster-cluster connections did not intersect voids we would 
have ended up with a larger overall percentage of filaments. 

Nevertheless, it is interesting to look at the abundance 
of filaments for a given separation between two clusters. Fig- 
ure]^ shows the fractions of cluster-cluster connections that 
contain a filament - straight or warped - as a function of 
the length of those connections. Very close cluster pairs are 
always connected by a filament, and about a third of all 
cluster pairs whose separation is between 15 /i^^ Mpc and 
20 h~^ Mpc are connected by a filament. 

The fact that very close pairs of clusters - with sep- 
arations of up to 5/i~^Mpc - are always connected by a 
filament can be explained by the presence of cluster infall 
regions. As discussed in, e.g., Diaferio & Geller 1997, the 
infall region of a cluster extends out to about three times 
its virial radius. What this means is that two very close 
clusters do not constitute two strictly separated systems as 
their infall regions overlap, and if they are gravitationally 
bound, the two clusters may eventually merge. Indeed, close 
cluster pairs show the presence of a filament between them 
(Dietrich et al. 20o|)- 

Figure |7] shows the length fractions of straight and 
warped filaments (strictly speaking, for warped filaments the 
length I is not the length of the filament but the length of 
the cluster-cluster connection). On the average, straight fil- 
aments are shorter than warped filaments. Two thirds of all 
straight filaments are in the 5 /i~^ Mpc to 15 h~^ Mpc range 
whereas warped filaments have a much broader distribution 
that extends out to fairly large separations. Clearly, tidal 
fields must play a role here. In fact the visual impression is 



that many warped filaments are bent towards another mass 
concentration. 

The most interesting properties of filaments are the 
amount of matter and its spatial distribution inside the fil- 
ament. We investigated these quantities for our sample of 
straight filaments as follows®. As noted above, all filaments 
were identified visually. In almost all cases the axes of the 
straight filaments deviate somewhat from the cluster-cluster 
axis. Therefore, for each filament we compute the actual axis 
by projecting all particles onto the plane perpendicular to 
the cluster-cluster axis and then computing the center of 
mass of the particles in that plane. The resulting center of 
mass is located where the filament's axis intersects with the 
plane. We then use this axis to compute the enclosed (radial) 
mass profile by averaging over all angles. 

The resulting profiles show a very interesting behaviour. 
There is a fairly large variation in the profiles close to the 
filaments' axes. However, for each filament we find that at 
some radius, rs, the profile starts following an r~^ power law 
closely. For each filament, we determine that radius, r^, by 
finding the part of the profile that can be fit by an r~^ power 
law. As an example. Figure |H| shows the radial overdensity 
profile of one filament. The bold vertical line marks r^. The 
existence of the radius Ts indicates that filaments have a 
well-defined edge. With the mass of a filament M/« (> r^) ~ 
const, but the volume growing as r^ (the length I = const.) 
the enclosed overdensity of the filament 5{< r) oc r~^. 

In Figure |5| we plot the distribution of the radii rs. The 
majority of filaments posess radii rs between 1.0 ft"'^ Mpc 
and 2.0 h~^ Mpc but there are also narrower and wider fila- 
ments. 

Given that filaments have well-defined edges how much 
matter is contained in a filament and does the amount of 
matter depend on the length of the filament? Figure 1101 
shows for each filament the enclosed overdensity at r^ versus 
the length of the filament. There is no obvious trend with 
length, but there is a fairly large scatter in the enclosed 
matter. There are even a few cases where the enclosed over- 
density exceeds 30. We visually inspected these cases and 
found that for every such case there is at least one very 
large halo present that lies on the filament's axis and whose 
mass is below the cluster threshold mass used for this work. 

The radial density profiles of filaments have quite in- 
teresting repercussions for observational efforts to find fila- 
ments. Given that filaments have well-defined edges with no 
dependences on length and with fairly large possible masses 
finding them observationally could be easier than previously 
thought. Either surveys using gravitational lensing or direct 
o bservations o f galaxies in the vicinity of galaxy clusters 
(lEbeling et al.li2004) appear to be very promising in this 
light. 

Figure UTI shows the averaged longitudinal overdensity 
profile of straight filaments. For these, we normalized the 
length of all filaments to unity and excluded cluster pairs 
with separations less than 5/i~^Mpc. We then computed 
the enclosed density of all material that is contained within 
2 h~^ Mpc from the filament's axis as a function of the posi- 
tion along the axis. The overdensity rises towards the clus- 



® Unfortunately, warped and irregular filaments cannot be inves- 
tigated in such a straightforward fashion. 
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Figure 12. Fractions of filaments for clusters: The plot shows the 
fractions of clusters that have 0, 1, 2, ... filaments (both straight 
and warped). Less than 10% of all clusters are not connected to 
a filament. 



Figure 14. For each cluster, this plot shows its number of fila- 
ments as a function of the distance to its fifth-nearest neighbour- 
ing cluster. 
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Figure 13. Average number of filaments for clusters as a function 
of cluster mass. Shown are clusters from the simulation (solid 
line) and clusters whose positions have been randomized (dotted 
line - see main text for details). Overplotted error bars show the 
standard deviation. The number of filaments clearly increases as 
clusters get more massive. 



ters and is constant at a value of around 7 in between the 
clusters. However, one needs to be careful with the profile for 
the following two reasons. First, individual filaments show 
quite a bit of irregular lumpiness when plotted this way. Av- 
eraging over the whole sample removes the lumpiness. Thus, 
while the averaged longitudinal profile is fairly smooth, indi- 
vidual examples look vastly different. And second, the pro- 
cedure we use neglects the fact that infall regions have dif- 
ferent sizes. We also plotted the averaged longitudinal pro- 
files of subsamples of filaments whose lengths are compa- 
rable. Because of the individual lumpiness of the filaments 
and because of the modest samples sizes the different sizes 
of the infall regions are lost because of the scatter of the 
data. Therefore, while averaging over all filaments longer 
than 5 h~^ Mpc is not ideal, our sample does not allow us to 
do more detailed studies. 



4-5.2 Clusters and filaments 

Since we started out with clusters to find filaments, it is quite 
natural to come back to them. The first obvious question to 
ask is how many filaments we can find per cluster. Figure 
1121 shows the distribution of filaments per cluster. The vast 
majority of clusters are connected to between one and four 
filaments. 

The average number of filaments per cluster is a some- 
what unsatisfactory quantity on its own. It gets more inter- 
esting when one looks at the mean number of filaments per 
cluster as a function of cluster mass - shown as the solid 
histogramme in Figure IT^ There is a clear trend with mass. 
More massive clusters are connected to more filaments. For 
the most massive clusters in the sample, the mean number 
of filaments is almost five. This picture agrees wit h the vi- 
sual im pre ssion from simulat i ons (s ee the images in I Jenkins! 
(1998) and lKauffmann et al.l il99iJ) ) where the backbone of 
large-scale structure is formed by the most massive objects. 

We tried to evaluate this result by randomizing the clus- 
ter positions and then looking for filaments again. The clus- 
ter positions were changed as follows. We first produced a 
smoothed version of the density field in the simulation vol- 
ume, using a Top Hat with a radius of 2 h~^ Mpc. We then 
assigned the clusters randomly to cells whose overdensity 
was five or more. This way, we made sure that the random- 
ized clusters ended up somewhere in the overdense parts of 
the simulation volume and not in a void. Having obtained 
this new set of clusters we re-did the investigation of the 
cluster-cluster connections using the same procedure as for 
the original cluster sample. The dotted histogramme in Fig- 
ure E] shows the result. As can be seen, the randomized 
clusters deviate from the original ones. The average numbers 
of filaments per cluster are lower than those of the original 
clusters and there is no trend with mass. 

Using a scatter plot. Figure IT^ shows the number of fil- 
aments per cluster versus the distance to the clusters' fifth- 
nearest neighbour. The latter is commonly used as a some- 
what crude measure of density. There is a trend for clusters 
with closer neighbours to have more filaments. This can be 
understood from the preceding: We have already seen that 
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Figure 15. Longitudinal velocity profile of straight filaments, 
averaged over straight filaments that are longer than 5h~^ Mpc 
(solid line) and five longitudinal velocity profiles of individual 
filaments (dotted lines). Shown is the enclosed mean velocity as 
a function of the positions along the cluster— cluster axis for all 
material that is contained with 2/i~^Mpc from the axis. The 
velocities are negative for material that moves towards the cluster 
on the lefthand side and positive for material that moves towards 
the cluster on the righthand side. 




Figure 16. Transversal velocity profiles of the five straight fila- 
ments, which were shown in Figure [T51 (dotted lines). Shown is 
the enclosed mean velocity as a function of the positions along 
the cluster— cluster axis for all material that is contained with 
2 h~^ Mpc from the axis. 



more massive clusters connect to more filaments. And we 
also know that in CDM universes more massive objects are 
more strongly clustered than less massive objects. 



4.5.3 Velocity fields 

As a consequence of the role of fila ments in the context o f 
the formation of galaxy clusters (see lColberg et all lll999l) ) 
it is most natural to assume that the velocity field at any 
given point inside a filament is determined by the distance 
to the two clusters. In particular, material probably moves 
towards the cluster it is nearest to. We use the straight fil- 
aments in our sample to investigate this. In order to avoid 



Figure 17. Fractions of clusters for which there is a filament 
in the connection to the nearest neighbour as a function of the 
distance. 



contamination through very close clusters pairs, which might 
be merging or whose infall regions overlap, we only consider 
cluster pairs that are separated by 5h~^ Mpc plus their re- 
spective virial radii. 

Figure \TEl shows the averaged longitudinal velocity pro- 
file of straight filaments (solid line) and the longitudinal 
velocity profiles of five individual straight filaments (dot- 
ted fines). For the profiles, we normalized the length of all 
filaments to unity. We then computed the enclosed mean 
velocity of all material that is contained within 2/i~^Mpc 
from the filament's axis as a function of the position along 
the axis. The velocities are negative for material that moves 
towards the cluster on the lefthand side and positive for ma- 
terial that moves towards the cluster on the righthand side. 
As can be seen, the average profile very clearly shows the 
gravitational domains of the two clusters. On the average, 
material tends to move towards the nearest cluster. The in- 
dividual profiles, however, show a fair amount of scatter. 
While the overall trend is the same as for the average pro- 
file, individual longitudinal velocity profiles are not nearly as 
smooth as the average. This fact can be understood from the 
clumpy structure of filaments. As seen above, the material in 
filaments is not distributed smoothly. Instead, single haloes 
determine the structure of a filament, with fairly large dif- 
ferences between individual filaments. Figure [TBI shows this 
very clearly. Here, we plot the transversal velocity profiles 
of the five straight filaments, which were shown in Figure 
1151 The scatter between the filaments is very large, and it is 
eas y to make out ind i vidua l haloes. 

lEisenstein et alJ (|l993) proposed a method to measure 
the mass of filaments by studying their transverse veloc- 
ity dispersion. We intend to investigate this method in an 
upcoming paper where we make use of a much larger simu- 
lation, which, with about 600 times as many particles in a 
volume about 44 the size of the simulation used here, has 
many more filaments and a much better mass resolution. 



4.5.4 Cluster-cluster alignments 

For a long time, there has been discussion on whether 
clusters of galaxies are aligned with their neighbours 
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IChambers et alJ ( 2000l) and iPlionis "fc Basilakod J2002D for 
recent updates). lOmiora fc Thomaa (|200fl l studied align- 
ments of clusters in large simulations and found "strongly 
significant alignments" for se parations of up to 30 h~^ Mpc 
in the ACDM model (also see lFaltenbacher et al]l2n02fl . As 
an explanation of align ment - if the effect exists - filaments 
have been brought up. I West et al.l (Il995[) sug gested clus- 
ter formation along filaments ( compare IColberg et al. 1999 ) 
had implications on the orientations of clusters. There are 
two possible explanations for this effect. First, the primor- 
dial density field pre-determines the directions from which 
matter falls into clusters. There is a positive correlation 
between the inertia tensor of a clust er and its surround- 
in g tidal field (see iBon d et alJ Il99ff). S econd, as shown 
in Ivan Haarlem fc van de Wevgaerti il993r . clusters tend to 
orient themselves towards the direction of l ast matter in- 
fall. T his finding, combined with the results of lColberg et alJ 
l|l99^ suggests that we can expect neighbouring clusters to 
be aligned if the formation of clusters along filaments is the 
dominant factor, which decides about alignment, and if there 
actually is a filament between the clusters. 

We have already seen earlier ("Section l4.5.2ll that not all 
neighbouring clusters are connected by a filament. Because 
alignment studies often focus on the nearest neighbour we 
want to shed additional light on this point in this context. 
Figure 1171 shows the fractions of clusters for which there is 
a filament in the connection to the nearest neighbour plot- 
ted against the length of the connection. Note that this is 
somewhat of a variation of Figure |5] Instead of looking at 
all cluster-cluster connections, we look only at those for the 
nearest neighbour of each cluster. As before, we see that very 
close pairs of clusters are always connected by a filament. For 
larger separations, the fraction drops, but even for separa- 
tions between 10 h~^ Mpc and 15 h~^ Mpc, the likelihood of 
finding a filament towards the nearest neighbour of a cluster 
is around 75%. Therefore, if the alignment of neighbouring 
clusters is caused by matter infall from filaments, there is a 
very high chance of finding cluster alignment up to separa- 
tions of 15 h~^ Mpc. Going further out, the likelihood drops 
quite steeply. 

There is another consequence of what we just discussed. 
If the formation of clusters along filaments is responsible for 
alignment, the non-detection of alignment of neighbouring 
clusters - separated by 5h~^ Mpc or more - does not nec- 
essarily mean that there is no such effect. It could simply 
mean that there is no filament between the clusters. 



4-5.5 Clusters and sheets 

As noted earlier, we were able to find a few sheet-like con- 
nections between clusters. We want to emphasize first that, 
as for the case of the filaments we discuss, we are only look- 
ing at sheet-like configurations of matter between neigh- 
bouring clusters. Whether or not some clusters lie in larger 
sheets - in galaxy redshift surveys these sheets are usually 
called walls - is beyond the scope of this work. We also 
note that some of the sheets could be mere projection ef- 
fects and we might have missed some sheets and classified 
them as projection effects. Given this uncertainty and given 
the small number of sheets we do not look at sheets in too 
much detail. 



Sheets appear to be much rarer than filaments. We 
found a total of less than two dozen sheets. In almost all 
cases there was only one sheet per cluster. The mean length 
of the connections that contain a sheet is 27.4 h~^ Mpc. We 
did not find very short sheets - with lengths smaller than 
10 h~^ Mpc - or very long sheets - with lengths exceeding 
45 /i"^ Mpc. 

Sheets appear to have lower surface densities than fil- 
aments. One might argue that the mass resolution of the 
simulation is not good enough to properly address the issue 
of how many sheets there really are. While we think that the 
mass resolution of our simulation definitely is high enough 
and that sheets indeed are not as common as filaments we 
want to note that we intend to re-address this issue in later 
work that will make use of a larger simulation with much 
higher mass resolution. 



5 SUMMARY 

Using a high-resolution N-body simulation of cosmic struc- 
ture formation in a ACDM Universe, we studied the mate- 
rial between pairs of clusters that correspond to R=0 Abell 
clusters. Our main results can be summarized as follows: 

• Whereas sheets appear to be fairly rare, filaments be- 
tween clusters are very common. The likelihood of finding a 
filament between two neighbouring clusters increases as the 
separation of the clusters decreases. Very close pairs of clus- 
ters that are separated by 5h~^ Mpc are always connected 
by a filament. Warped or irregular filaments are more com- 
mon than straight filaments - due to the presence of tidal 
fields between clusters. The longer the cluster-cluster con- 
nection the higher the likelihood is of finding a warped or 
irregular filament rather than a stright one. 

• We investigated straight filaments in more detail to look 
at the amount of matter in those filaments. Filaments are 
very lumpy objects. The longitudinal overdensity profile - 
measured along the axis that connects the two clusters - 
clearly shows this. Towards the two clusters the overdensity 
rises as the cluster infall regions are reached. The radial 
enclosed overdensity profiles of filaments show a well-defined 
radius at which the profiles follows an r~^ profile. For the 
majority of filaments, this radius is between 1.0 h~^ Mpc and 
2.0 h~^ Mpc but there are also narrower and wider filaments. 
The enclosed overdensity inside this radius varies between a 
few times up to 25 times mean density or more. All high- 
density cases could be visually identified as containing large 
haloes. 

The results from the density profiles indicate that find- 
ing filaments observationally might be somewhat easier than 
previously thought especially if the line of sight is aligned 
with the axis of a filament. Finding filaments that are per- 
pendicular to the line of sight is trickier because of the lumpi- 
ness of filaments. It is probably most promising to look at the 
immediate vicinities of clusters where filaments have higher 
densities. 

• The majority of all clusters possess between one and 
four filaments. There is a very clear trend with mass. More 
massive clusters on average have more filaments. This result 
supports the general view that the most massive clusters sit 
at the intersections of the backbone of large-scale structure. 
There also is a weak trend for clusters in denser regions to 
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have more filaments. The last statement is not completely 
independet of the former one. More massive clusters can be 
found in denser regions and are also more clustered. 

• The velocity field in a filament is dominated by the two 
clusters. While there is a large scatter between individual 
filaments, material tends to move towards the clusters that 
is closest. 

• Filaments have been used to explain alignments be- 
tween neighbouring clusters. We find that there is a very 
high likelihood of finding a filament in the connection to 
nearest neighbour of a cluster for separations of up to 
15 h~^ Mpc. If filaments are indeed responsible for align- 
ments one would expect to find closer pairs of clusters to 
be aligned. However, even at separations of between 5 and 
15 h'^ Mpc there are pairs of clusters that are not connected 
by a filament. Close clusters that are not alligned thus could 
indicate that there is no filament between them. 

• We find that the fraction of matter configurations that 
appear to be filamentary due to projection effects is about 
the same as the fraction of genuine filaments. This effect 
has to be taken into account when investigating the dis- 
tribution of matter between neighbouring clusters observa- 
tionally. The contribution of projection effects become non- 
negligible at separations of 10 h~^ Mpc. 

The result s of t his paper have found application in 
IPimbblet et alJ (|200J), who analyze the frequency and dis- 
tribution of intercluster galaxy filaments selected from the 
2dF Galaxy Redshift Survey. 
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